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Abstract. We discuss universal properties of higher order cumulants of net baryon 
number fluctuations and point out their relevance for the analysis of freeze-out and 
critical conditions in heavy ion collisions at LHC and RHIC. 



1. Chemical freeze-out, the hadron resonance gas and lattice QCD 



Higher order moments, or more accurately higher order cumulants, of conserved charges 
are central observables analyzed in the ongoing low energy runs at RHIC. It also is 



CO . quite straightforward to analyze their thermal properties in equilibrium thermodynamics 

00 ■ 

of QCD, e.g. by performing lattice QCD calculations. Ratios of cumulants of 
different order are particularly well suited for a comparison with experiment as they 
are independent of the interaction volume. Lattice QCD calculations have shown that 
such ratios, e.g. ratios of baryon number, electric charge or strangeness fluctuations, 
are quite sensitive probes for detecting critical behavior in QCD. They are sensitive 
to universal scaling properties at vanishing as well as non-vanishing baryon chemical 
potential (/is) and directly reflect the internal degrees of freedom that are carriers of 
the corresponding conserved charge [TJ [2] in a thermal medium. These ratios change 
rapidly in the crossover region corresponding to the chiral transition in QCD and reflect 
the change from hadronic to partonic degrees of freedom [3]. 

In the low temperature phase of QCD ratios of cumulants seem to be well described 
by a hadron resonance gas model (HRG) [U El H]. The HRG is also very successful in 
describing the thermal conditions which characterize the chemical freeze-out of hadron 
species. The freeze-out temperature and its dependence on /i^ is found to be close to 
the pseudo-critical temperature T pc for the QCD transition. We recently pointed out 
that also ratios of cumulants of net baryon number fluctuations, as measured by STAR 
at different beam energies [5], agree well with HRG model calculations on the chemical 
freeze-out curve j6]. This is shown in Fig. [TJ At the same time these results are also 
consistent with lattice QCD calculations when temperature and chemical potential in 
these calculations is chosen to agree with the freeze-out parameters [TJ [8] . This also is 
shown in Fig. [TJ 
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Figure 1. Ratios of cumulants of net baryon number fluctuations measured by 
STAR [5] and calculated in a HRG model [6]. The lattice QCD results [8] have been 
determined from a next-to-leading order Taylor expansion of cumulants calculated at 
the values of the freeze-out chemical potential determined from the HRG model. The 
freeze-out temperature at fis = has been assumed to coincide with the crossover 
temperature determined in lattice calculations. 



The HRG model, of course, is not sensitive to critical behavior. The agreement of 
current lattice calculations with HRG model calculations of ratios of cumulants including 
up to fourth order fluctuations thus also indicates that they present so far no compelling 
evidence for critical behavior. We recently argued that a determination of sixth order 
cumulants may be particularly helpful in this respect as they clearly deviate from HRG 
model calculations in the crossover region of QCD and start to become sensitive to 
critical behavior earlier in the hadronic phase [9] . In the following we will point out some 
generic features of higher order cumulants in QCD that arise from universal properties 
of QCD close to the chiral limit of vanishing light quark masses, m q = m u = nid = 0. 

2. 0(4) universality 

At vanishing baryon chemical potential (/ie), as well as in a certain range < [ib < A^j 
QCD is expected to undergo a second order phase transition for m q = 0. In this 
entire range of \xb values the chiral phase transition belongs to the universality class 
of 3-dimensional, 0(4) symmetric spin models. In the vicinity of the chiral phase 
transition temperature (T c (/2b)) thermodynamic quantities show universal properties 
that are controlled by the singular part, //, of the free energy. For a fixed value of the 
chemical potential JIb < A 4 ^, close to T c (Jj,b) the free energy may be parametrized as, 



with z = t/h l ^ s denoting the particular scaling combination of the reduced temperature 
t ~ (T - T c (J1 b ))/T c {JXb) + k b {{(iB/T) 2 - (/2b/T) 2 ) and the symmetry breaking 
parameter h, which for QCD is taken to be the ratio of light and strange quark 
masses, h ~ m q /m s ; a, (3, S denote critical exponents of the 3-d, 0(4) universality 
class. Cumulants of net baryon number fluctuations are then obtained as derivatives of 
f{T,fiB,rn q ) with respect to Jib = Hb/T, i.e., Xn ~ d n f/dfi B - Close to the chiral limit 



f(T,fi B ,m q ) = h 2 



ff(z) + fr(T,H B ,m q ) 
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Figure 2. The third derivative of the singular part of the free energy in theories 
belonging to the 3-d, 0(4) universality class (left) and its contribution to third or higher 
order cumulants (see Eq. 2) (right), ho and zq are non- universal scale parameters. 

and for n > 3 these derivatives are dominated by the singular part, 

f -mf- a ~ n/2)/ps ff /2 \z) , for fx B /T = 0, and n even 

X " ~ ( - Jmf-^ffiz) , for Hb/T > ^ 

The first derivative of the scaling function that leads in the chiral limit to a divergent 
cumulant at T c (/ib) and, in fact, leads to a change of sign of \n a ^ T c {Jib) is obtained 
from the third derivative of the scaling function fi 3 \ In Fig. [2] (left) we show this scaling 
function for the 3-d, 0(4) universality class which recently has been extracted from high 
statistics spin model calculations [10] . The right hand part of this figure shows how the 
singular contribution strengthens as the symmetry breaking parameter h is reduced. 

It is evident from Fig. fright) that the universal scaling properties of cumulants 
will, for sufficiently small values of the quark mass, induce negative values for cumulants 
already for n = 3. However, as shown in Eq. El the contribution from the singular part is 
weighted by a factor proportional to (//b/T) 3 . Its contribution to the total value of the 
cumulant may thus be small for small values of Hb/T and non- zero values of the quark 
mass. Nonetheless, negative values for 3 rd and 4 th order moments have been found in 
model calculations fTTjljl 

The first cumulant, which is not parametrically suppressed by powers of hb/T and 
thus stays finite even at /i# = is the sixth order cumulant of net baryon number 
fluctuations. This cumulant changes sign in the crossover region of the QCD transition. 
Although the temperature at which the sixth order cumulant changes sign is not a 
universal quantity there are strong indications from lattice jl] as well as model [9] 
calculations that the change of sign occurs close to the chiral transition temperature. 
This is shown in Fig. [3j In fact, for sufficiently small values of the quark mass, i.e., 
in the 0(4) scaling regime, the location of the minimum of the sixth order cumulant 

X Close to the chiral critical point the negative values of fourth order cumulants that arise from 0(4) 
criticality compete with similar effects that arise from Z(2) critical behavior and also lead to negative 
fourth order cumulants [12) . 
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Figure 3. Sixth order cumulants of net baryon number, electric charge and strangeness 
fluctuations calculated in (2+l)-flavor QCD at ^ B = gj [5] (left) and in the PQM 
model (right) [S]. 



and its position relative to the pseudo-critical temperature for the chiral transition is 

(3) 

controlled by the location of the maximum of n and its location relative to the peak 
in the chiral susceptibility. The latter appears at a somewhat higher temperature [TO]. 
For [1b/T > the onset of negative values for xi indeed follows the crossover line for 
the QCD transition P EL2|. 

3. Conclusions 

Sixth order cumulants are thus expected to change sign at a temperature below the 
(pseudo-critical) chiral transition temperature. This effect should become visible even 
at the LHC and the highest RHIC energy if chemical freeze out indeed occurs close to the 
QCD transition temperature and if higher moments probe these freeze out conditions. 
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